Direct experimental evidence of free fermion antibunching 
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Fermion antibunching was observed on a beam of free noninteracting neutrons. A monochromatic 
beam of thermal neutrons was first split by a graphite single crystal, then fed to two detectors, 
displaying a reduced coincidence rate. The result is a fermionic complement to the Hanbury Brown 
and Twiss effect for photons. 
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Over the past three decades, the research on the foun- 
dations of quantum mechanics has been enriched by 
many experiments on thermal neutrons, in particular by 
several enlightening results about the coherence prop- 
erties and the physical nature of the wave function de- 
scribing the behavior of a massive particle [1] . A general 
property of fermions is that of being characterized by an 
antisymmetric wave function: the second-order correla- 
tion function of a fermion gas exhibits an anticorrelation 
in the intensity fluctuations, in particular interference in 
the coincidence distributions of identical particles. The 
present Letter describes a new contribution in this field: 
an experiment on thermal neutrons that brings to light 
the fermion antibunching effect in a beam of free nonin- 
teracting particles. A monochromatic beam of thermal 
neutrons was first split by a graphite single crystal, then 
fed to two detectors, displaying a reduced coincidence 
rate. The result is a fermionic complement to the seminal 
Hanbury Brown and Twiss effect for bosons (photons) [2]. 

The consequences of antisymmetry are well known in 
condensed matter physics, where the electronic states dis- 
play a strong quantum entanglement and are confined 
within the Fermi surface. Interesting experiments with 
electron beams have confirmed these effects [3-5]. In 
the case of almost free particles, an anticorrelation was 
observed in the coincidence spectrum of neutrons from 
compound-nuclear reaction at small relative momentum 
[6, 7]. However, such a physical system is not a good 
representative sample of a statistical ensemble of non- 
interacting identical fermions. A monochromatic beam of 
thermal neutrons from a nuclear reactor represents much 
better a statistical ensemble of free particles. Neverthe- 
less, the observation of thermal-neutron antibunching by 
means of coincidence measurements on such beams with 
the available instrumentation did not appear to be fea- 
sible up to now, mainly because the mean number of 
fermions obtainable per unit cell of phase-space, to which 
the signal-to-noise ratio is proportional, was so low that 
a measurement time of several years was estimated [8] . 



We shall show below that, with present-day available 
advanced instrumentation, a very accurately designed 
setup and a precise knowledge of the statistical proper- 
ties of the neutron source, the experiment is feasible. In 
this Letter we shall describe some measurements carried 
out at the Institute Laue Langevin, Grenoble, France. 

How can one directly bring to light an anticorrelation 
effect in a neutron beam? In a gas of fermions there 
is a certain tendency for particles of the same spin to 
avoid each other, a tendency arising from the exchange 
antisymmetry of the wave function: two fermions in the 
same spin state cannot occupy at the same time the same 
point in space, and therefore the probability amplitude 
for their being close together must be small. We just want 
to observe such an effect in a beam of thermal neutrons. 

Let us start by considering the conceptual scheme of 
our experiment, which is schematically represented in 
Fig. 1, and is a massive particle analogue of the seminal 
optical Hanbury-Brown and Twiss experiment [2] , which 
yielded the first direct observation of the bunching effect 
in light beams and is a direct consequence of the symmet- 
ric wave function of a bosonic state. The semi-classical 
and quantum interpretations of this experiment are very 
clearly discussed in classic textbooks [9, 10]. 

In the present experiment, a monochromatic beam of 
thermal neutrons is split into two components by a py- 
rolytic graphite single crystal that produces a transmit- 
ted and a reflected beam. The intensity of each com- 
ponent is measured by a detector, and the coincidence 
rate of the outcomes at the two detectors is recorded 
by the coincidence counter as a function of their rela- 
tive distance from the splitting crystal. Of course, also 
accidental coincidences occur in the apparatus and their 
rate must be subtracted from the total rate. The av- 
erage value of these rates over a long enough period of 
observation provides the quantity to be analyzed. 

The nature of the emission of thermal neutrons in the 
source is expected to be Poissonian, so that there is a 
small but finite probability of having two neutrons within 
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FIG. 1: Schematic drawing of the experimental setup: 
M=monochromator, S=beam splitter, Di l 2=detectors, 
C=coincidence counter, DAE=Data Acquisition Electronics. 
The two detectors can be positioned at the same distance 
from S, and one of them can be moved across this distance. 
The collimators are not shown. 



the detection time interval of the apparatus. With ref- 
erence to Fig. 1, for an average total rate n of neutrons 
impinging on the splitter crystal S and a total measur- 
ing time To, the predicted number of coincidences N c 
measured at the two detectors can be calculated from 
the joint probability that a neutron is transmitted while 
a second neutron is reflected from S. By neglecting the 
fermion antibunching, one readily obtains: 
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where n t and n d — n — n t are the average rates of 
the transmitted and reflected beam respectively, that 
have been assumed to be constant during the acquisi- 
tion time To, and we assumed a short enough coincidence 
window t w /2, such that nr w <C 1 (in our experiment 
nr w ~ 10~ 2 -=- 1CP 3 ). When the distances of the two 
detectors from the beam splitter, Sui and Sd2, are dif- 
ferent enough, Eq. (1) gives the expected number of ran- 
dom coincidences, since the measured coincidences are 
associated only with the simultaneous detection of two 
particles emerging from the splitter at different times. 

Let us now qualitatively analyze the expected conse- 
quences of antibunching. Again referring to an ideal ex- 
periment, let Sdi = <Sd2 and r c be the coherence time of 
the neutron wave packet. It should be remarked that the 
coherence time of the wave packet is defined by its en- 
ergy distribution and it is longer when the energy width is 
small. Equation (1) can be modified to yield the expected 
number of correlated coincidences: observe that these are 
only those due to two neutrons with different spins that 
emerge from the beam splitter at the same time, because 
two neutrons with the same spin, due to the Pauli ex- 



clusion principle, cannot impinge on the beam splitter at 
the same time. If the incident beam is spin-unpolarized, 
it is equally likely that a neutron pair will either occur in 
one of the three triplet states or in the singlet state, i.e. 
the triplet states will occur 3/4 of the time. Thus, the av- 
erage number of coincidences expected for Sui = Sd2, as 
a consequence of fermion antisymmetry, is reduced from 
that of Eq. (1) by the following quantity 

-AA fa = -^T c n t n d T , (2) 

where -1/2 = -3/4 + 1/4, -3/4 being due to the an- 
tibunching of the triplet states and 1/4 to the bunching 
of the singlet state. Such a depression of the coincidence 
rate involves a two-particle state, and it is essential that 
both members of the pair be detected. Since the rela- 
tive directions of the two particles may not be known 
in advance, one might think that something near a 47r 
detector might generally be needed. But in the present 
experiment the two-particle state to be tested is only that 
emerging (within a small solid angle) from the collimator 
and the monochromator, and the expression r c n t rid in 
Eq. (2), is simply the rate of such emerging state. This 
must be taken into account if one plans to perform a 
neutron-spin test of the Bell inequality [11]. 

Of course, in a real coincidence experiment, one must 
take t w much longer than r c , in order to account for 
various instrumental effects that force one to broaden 
the coincidence window. Actually, two particles arriving 
at the beam splitter at the same time may be absorbed 
and recorded at the two detectors within a rather long 
time interval td, because of the finite thickness of the 
beam splitter and detectors and of small differences in 
their speeds. Moreover, the finite detection resolution 
td < T w tends to reduce the value of Nc^ in Eq. (1) by 
a factor t w /td, namely 
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making the use of intrinsically fast detectors highly desir- 
able. Clearly, as far as td ^> t c , AN[ a in Eq. (2) remains 
unchanged. The experimental data require therefore a 
careful analysis, as we shall discuss in detail in the fol- 
lowing. 

In order to detect the expected fermion antibunching 
effect, we have performed an optimized experiment based 
on the general scheme of Fig. 1. The main limitations of 
the experiment arise from the random fluctuations, which 
can mask the difference signal of Eq. (2). We assume that 
there are no accidental coincidences due to non-random 
processes. Therefore the expected random fluctuations 
are those due to the intrinsic statistics of the number of 
measured coincidences. The expected root mean square 
(rms) fluctuation of the total number of coincidences is 
given by: 



AA C = \J T D n t n d T . 



(4) 
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In order to detect a signal it is necessary that the noise 
to signal ratio, 
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be smaller than unity. Looking at Eq. (5) we see that the 
noise to signal ratio decreases when the coherence time 
of the incoming beam is long, so that the experiment 
should be performed by employing the most monochro- 
matic available beam. We chose to use the IN10 beam 
line [12] , which produces a monochromatic beam by using 
an almost perfect Si(lll) single crystal in the backscat- 
tcring configuration. This monochromator produces a 
flux n ~ 2000 sec -1 , at an energy E ~ 2.08 meV with an 
energy spread AE < 0.13 ^eV (FWHM= 0.3 ficV). The 
coherence time of the incoming neutron beam is there- 
fore r c > H/2AE ~ 2.5 ns. Considering that the neu- 
tron speed v is about 630 m/s, the neutron coherence 
length is larger than 1.5 /mi, a very small value. It is 
clear that both the beam splitter and the detectors must 
be as thin as possible, in order to reduce any additional 
spread of the signal. We therefore employed a 0.3 mm 
thick graphite crystal as beam splitter. At the wave- 
length of the present experiment the crystal has a good 
reflectivity so that the transmitted and diffracted beam 
are of the same order of magnitude. It should also be re- 
marked that the possible velocity difference between the 
two particles of a pair, originating from the energy spread 
of the monochromator, contributes a negligible difference 
in the corresponding time of flight along the 40 cm path 
from the splitter to the detectors. 

Two different detection systems were employed. The 
first one was based on two squashed 3 He 2 mm thick, 
1.2 cm wide and 10 cm high detectors. The second one 
was based on two scintillator detectors having thickness 
0.2 mm, width 1.5 cm and height 5 cm. The scintillator 
was a 6 Li 98% enriched ZnS glass, directly coupled to a 
5 cm diameter fast photomultiplier. The shaping time 
was about 2 fis in the case of 3 He detectors, while it was 
0.3 /is in the case of the scintillators. The coincidences 
were measured within a time window t w /2 of ±10 /us in 
the case of gas detectors and ±0.8 /is for the scintillators. 

Using this arrangement we have been able to perform 
two meaningful determinations of the actual antibunch- 
ing effect in the incoming neutron beam. In order to 
do this, one detector was kept at a fixed position, at a 
distance of 40 cm from the graphite splitter along the 
diffracted beam, while the other detector was scanned 
through the transmitted beam, at approximately the 
same distance. For gas detectors, whose spatial reso- 
lution is of order 2 mm, we used a coarse translation step 
of 1 mm, while for scintillator detectors the translation 
step was 0.2 mm. 

The data acquisition took several days; we therefore 
had to take into account the effect of the incoming beam 
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FIG. 2: Number of coincidences N c as a function of the path 
difference x = Sui — S02 of the detectors from the splitter. 
Above: scintillator detectors, translation step 0.2 mm; below: 
gas detector, translation step 1mm. The dip appears at xo ^ 
2 mm due to calibration. The parameters td,t c , iVi 00 ' and 
xo, by substituting t = (x — xo)/v in Eq. (9), are determined 
by a best fit (see text). 



fluctuations. As can be seen from Eq. (3), such an ef- 
fect is nonlinear and directly related to the instantaneous 
value of the incoming beam rate. Since the instantaneous 
rate cannot be measured with adequate accuracy, one can 
perform a correction of the actual data by assuming that 
the beam fluctuations are small. In such a case, assuming 
that the acquisition time T is much longer and the detec- 
tion window t w much shorter than the correlation time 
of the noise and neglecting second order effects, Eqs. (2) 
and (3) yield 
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where N c is the actual number of coincidences, N c the 
number of coincidences detected in an (ideal) experiment 
with a constant rate on both the transmitted and incom- 
ing beam, and Nt and Nd are the actual numbers of neu- 
trons detected on the two beams. All the quantities in 
Eq. (6) are functions of time, but vary on timescales that 
are much slower than the time window t w . 
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The experimental data collected with the two detect- 
ing systems were corrected according to Eq. (6) and the 
results are reported in Fig. 2. In both cases a small dip 
is observed in the number of coincidences detected as 
a function of the relative distance of the two detectors 
from the beam splitter. We attribute this small dip to 
the antibunching effect due to the Fermion nature of the 
neutron. 

It is interesting to perform a quantitative analysis, us- 
ing the experimentally observed width of the dip in order 
to get an estimate of the coherence time of the incoming 
beam neutron wave packet. As a spinoff, this will yield 
a consistency check of our experimental results. Let us 
first consider the global response function of our exper- 
imental arrangement. We assume that the neutrons are 
diffracted within the thickness of the beam splitter ac- 
cording to the secondary extinction law [13] and are ab- 
sorbed by the detectors, within a few mean free paths 
inside the absorbing medium. For simplicity, let us as- 
sume that the total response function of the apparatus 
be (a normalized) Gaussian 

with a characteristic response time td S> t c , where t c 
is the coherence time defined in Eq. (2). The global re- 
sponse time td of the apparatus will be obtained by fit- 
ting the experimental data: it is however expected to be 
close to the shaping time of the detectors (2 (is for 3 Hc 
detectors, 0.3 (is for scintillators). 

The neutron pair correlation function, describing the 
antibunching effect, will also be taken to be Gaussian 

C(i) = l-iexp(-^), (8) 

the factor 1/2=3/4-1/4 being due to the difference 
between the triplet and the singlet contributions [see 
Eq. (2)]. The total number of counts is therefore given 
by the convolution of these two functions 

P| = [fl.C](<) = I dsR(,)C(t-,) 

This must be compared with the observed number of co- 
incidences in Fig. 2. Looking at the experimental data, 
we see that in both cases there is a small but apprecia- 
ble dip, which is broader in the case of the experiment 
performed using the (thicker) gas detectors, as expected. 
The above formula implies that the width of the dip is 
td, its depth being t c /2td = AN^/Nc 00 ^, in agreement 
with Eqs. (2) and (3). 

An accurate fit yields td = 1.3 ± 0.4 (is (with x = 
2.1 ± 0.2 mm and A c (oo) = 34720 ± 44) for 3 He detectors, 
and T B = 0.33 ± 0.07 (is (with x = 1.93 ±0.02 mm and 



N (°°> = 993.7 ± 0.6) for the scintillators. Moreover, 

t c = 20 ± 7 ns for 3 Hc detectors, (10) 
t c — 19 ± 3 ns for scintillators, (11) 

both values being fully consistent with each other and 
with the bound obtained by the energy spread of the 
beam (> 2.5 ns). The fitting curve is shown as a full line 
and is in very good agreement with the data. Note that 
the value td — 0.33 (is is in full accord with the nominal 
shaping time of the scintillator (0.3 (is), while the value 
td — 1.3 /xs is smaller than the nominal shaping time 
of the 3 He detectors (2 (is): this can be understood by 
remarking that 3 He detectors tend to absorb neutrons 
in the initial section. Finally, notice that one obtains 
results that are consistent with the above ones also by 
performing a convolution with more realistic (nongaus- 
sian) shape functions describing the response function of 
the detectors and of the beam splitter. 

It is useful to clarify in what sense this experiment 
performed with neutrons is complementary to its photon 
[2, 14] and electron [3-5] counterparts. The first, obvi- 
ous observation is that neutrons are fermions that are 
not affected by Coulomb interaction, that plays, by con- 
trast, an important role in condensed matter systems. 
Second, neutrons have very low phase-space densities, so 
that all the effects we have brought to light are due to 
two-particle correlations, three or more particle effects 
being completely negligible. This experiment, providing 
a firm experimental evidence for the Pauli exclusion prin- 
ciple, displaying its effects on free neutrons in real space, 
has a very basic importance, because it is directly related 
to the quantum mechanics of identical particles. 
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